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Abstract 

We derive a formula which relates the QCD potentials in momentum space and in 
position space in terms of the (5 function of the renormalization-group equation for the 
potential. This formula is used to study the theoretical uncertainties in the potential 
and in particular in its application to the determination of the pole mass nib when we 
use perturbative expansions. We demonstrate the existence of these uncertainties for the 
Richardson potential explicitly and then discuss the limited theoretical accuracy in the 
perturbative QCD potential. We conclude that a theoretical uncertainty of nib much 
below 100 MeV would not be achievable within perturbative QCD. 
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In this article we discuss a relation between the static QCD potential |IJ in momentum 
space 

V(,) = - Cf «, (!) 
where for quarks Cf = 4/3, and the corresponding potential in position spaceQ 

Eqs. (P and (||) are understood as the definitions of the couplings a v (q) and a v {l/r). In 
addition, they express the fact that these two functions are related because V(r) is the Fourier 
transform of V(q). We can obtain a useful expression for this functional relation 

a v (l/r) = F[a v (q)} (3) 

using the renormalizat ion-group equation of the coupling a v (q), and we will study the conse- 
quences of this relation. 

First we relate the coupling at a general scale, a v (q) , to its value at a specific renormalization 
scale \i. Let us express ot v (q) as a power series 

oo 

a v (q) = $(a v (//),t) = ^ c n (/i)t n where t = ln(> 2 / q 2 ) . (4) 

n=0 

It follows that 

a v(^) = $ ( a v(^)> 0) = cq(m)- ( 5 ) 

For static heavy quarks, the potential energy is a physical quantity. As a result, a v (q) obeys 
the renormalization-group equation 

2 d , , d<5>(a v ,t) d$(a v ,t) 

" av {q) = — ar~ + K } —^r = °' (6) 

where the j3 function for a v is defined by 

P V M = /i = -^XX„ [—) ■ (?) 

Since the coefficients c n (/i) in eq. (f|) can be expressed as partial derivatives of $ with respect 
to t, 



1 <9 n 

one can show using the renormalization-group equation 



(8) 

t=o 



1 ; 9 i 1 ; 9 

Cn(^) = ^| ( -Pv( a v) q^-J H a v,t= °) = —\ \~Pv( a v">]fa~) a v{v)- (9) 



1 We use the notation q — |q| and r = |r| 
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Combining this expression with eqs. (P and (H), we find the relation between the couplings 
a v (q) and a v (fi): 



a v (q) = exp 



-t(3 v (a v ) 







da. 



00 t n ( 

n=0 U - \ 







da. 



in+l 



From this expression, one finds 

cM = 0([a v (v)y 

Now we can perform the Fourier transform in eq. (|2|) using the following formula: 

d 3 q n e iq r d n 



J 



(2vr) 3 q 2 du 
where for — 1 < u < 1/2, 

^(r,/i, M )^/i 2 "| 



u=0 



d 3 q e iqr {fir'] 



2d 



with 



/(«) 



/tan(7ru) 

V ITU 

00 

£/»«" 
n=0 



exp 



(2vr) 3 (q 2 Y+ u 4vrr 



' ™ (2w) 2 " +1 



/(«) 



(10) 



;i2) 



(13) 



^ 2n + 1 



x + vr 2 ^ 2 + 8C(3)^ 3 + IQtt 4 ^ + / 4tt 2 C(3) + 32C(5) ' 
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Here. 



360 



r' — r e lE . 
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u° + 



(14) 



(15) 



7 E = 0.57721... is Euler's constant, and ( denotes the Riemann ( function.^ Higher order terms 
in the expansion of f(u) at u = can be easily obtained using e.g. Mathematica f|. It follows 
from eqs. (|T2| ) and ( pT3[ ) that 



F[t n ] 



Qn 

du n 



[e uy f(u)} 



= £ 

=0 k=0 



n 
k 



y 



k d n ~ k f(u) 



du 



n—k 



X^ n - f k 

/ j 77 J n—k V j 

k=0 K - 



u=0 



(16) 



2 Different expressions for ,F(r, fx, u) are derived in ||. One of those, 

rrruul _ (H 2 " r(| + u)r(|-u) 

j_ 4vr 2 r r(l + 2u) ; 
can be used to derive eqs. (Oh and (Il4|) via the formulas B 



va + z)va-z) = 



and 



In r(l + z) = I In " 1 



C0S(7TZ) 



00 2n+l 



2 sin(7rz) 
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where 



y = ln(/iV a ). (17) 

Then from eq. (|lO|) we obtain 

oo 771 Tj.nl / f) \ U 

a v (l/r) = F[a v (q)} = £ ~Y I -/?>v) «v(aO- (18) 

n=0 n - \ aQ V / 

In particular, for y — 0, which implies /x = 1/r', 

F[« n ](y = 0) = n!/ B . (19) 



Hence 



■ 

/ Q 



a 



V 



= £ /» "M«v) 7T— a v (g = 1/r'). (20) 



n=0 \ 9a v 



This equation is the specific representation of the functional relation (^) which we set out to 
derive. The function (3 v (a v ) is defined in eq. (^) and the coefficients /„ in eq. fll4|). 

Before studying consequences of eq. (p0|) for the perturbative QCD potential, let us discuss 
first a model which we can solve numerically and which describes the energy levels of quarkonia 
very well. The model that we consider is the Richardson potential || which in momentum 
space reads: 

Vniq) = -AttCf v = - ■ 21 

q 2 p q 2 \n(l + q 2 /Af l ) 

A very good description of the charmonium and bottomonium states is obtained using this 
potential for A R = 0.4 GeV and 

A) = 11 - \n f , (22) 

where nj = 3 is the number of the light quarks. In Fig.|l| the function a^ R \q) is plotted as 
the dash-dotted line for q between 1 and 10 GeV. The next-to-next-to-leading order (three- 
loop renormalization-group improved) QCD coupling a^\q) || for n/ = 4 active flavours and 
a TTs( m z) = 0.119 is also shown as the solid line.f] It is seen that in the range of q relevant for 
the position of T(IS') the two curves are fairly close. The QCD potential is more attractive 
which implies that the value of m& extracted from the mass of T(IS') is slightly larger for 
the next-to-next-to-leading order QCD than for the Richardson potential. In fact in a recent 
article || = 4.96 GeV is obtained^ with an uncertainty of about 100 MeV arising from the 
uncertainty in the input value of ajjg{mz)- For the Richardson potential mj = 4.88 GeV is 
obtained. All these attractive features suggest that at scales of a few GeV the model closely 
resembles the true QCD potential. 

3 We used the three-loop renormalization-group equation to evolve the MS coupling, and Eq. (9) of B 
without the a 3 -term to match the 4-flavour to the 5-flavour theory at the matching scale 5 GeV, and m& (pole) = 
4.88 GeV. The result is not very sensitive to varying this matching scale. 

4 This value is different from what B gives as a final value of mj,. It corresponds to the result from the 
two-loop static QCD potential, which does not include a shift in mj, of about 50 MeV due to relativistic and 
leading non-perturbative corrections. 
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Figure 1: Comparison of the momentum-space couplings: aty(q) - solid, a\y'(q) - dash-dotted, 
and a n (q) - dashed line. 



For this model, we study the validity of the relation ( p0|) when it is expressed as a pertur- 
bative series of a v in the asymptotic region where the coupling is small. In order to define the 
perturbative coupling of this model unambiguously in both momentum and position spaces, we 
subtract the confining part of the Richardson potential 



v R ( q ) = v con{ (q)-^c F ^-. 

q z 



(23) 



The confining part Konf(<?) corresponds to a linear confining potential in position space 

VcoAr) = ^A|r. (24) 

The coupling 



An 
To 



ln(l + g7A|) q* 



(25) 



is considered as a perturbative coupling in momentum space. The coupling a H (q) is plotted as 
the dashed line in Fig.|l]. The corresponding coupling in position space is given by the following 
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0.04958339 
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0.04257994 
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Table 1: Comparison of the coupling a R (l/r) for the Richardson potential and the best ap- 
proximation obtained by truncating the asymptotic series ( p9|) at n = N. 



expression ||: 

r°° du e _ARrM 



. . . 2ir 



1-4 



1 U ln 2 (M 2 - 1) + 7T 2 



(26) 



which can be calculated numerically to a very high precision. In particular the accuracy of our 
numerical calculations was 10" 10 . From eq. (|25| ) we calculate the function 



( a ) = q 2 ^ (g) = - + & [ ^ ~ " r]2 (o 7) 
q dq* 4tA + 47t 1 + g 2 /A^' 1 J 



where r/ 2 /A 2 l is understood as a function of a R . For g ^> Ar the second term in the expression 
for /3r is negligible. Then our problem becomes very simple: there is only one term in the sum 
in eq. (^) corresponding to 

Pr,o = Po- (28) 
Instead of eq. (^0]) we have a much simpler relation 

a R (l/r) ~ f>J 9 ) (^^) n\ f n , (29) 



n=0 



with q = 1/r'. The symbol "~" indicates that in our derivation non-perturbative higher twist 
contributions to /3r (terms suppressed by powers of e~ 47T ^ /3 ° aR ^) have been neglected. Strictly 
speaking, the relation Q2"9"D may be valid only for q ^> A R i.e. for small a R . We have computed 
the first twelve coefficients f n and checked that all are positive. For a given 1/r', the terms 
in fl2"§D decrease with increasing n for n < n min and increase for n > n min . The value n m i n 
corresponding to the minimal contribution increases with growing 1/r'. The best approximation 
is obtained when the series in (|29|) is truncated at AT ~ n mm . In Table the values are given 
of the functions a-n(l/r) computed numerically from eq. (|26|) and a^(l/r) obtained from the 



5 



asymptotic expansion fl2"9|) truncated at n = N. The coefficient (3q is evaluated for rif = 3 and 
1/r is varied between 10 and 10 7 GeV. It is evident that for large values of 1/r a very good 
approximation is obtained using the truncated series a-R^l/r). However, at 1/r = 10 GeV, 
^min = 1 and the quality of approximation practically does not improve when instead of 

a R (l/r)~a R (g=l/r') (30) 

the formula including the 0(a^) term is used. At 1/r = 20 GeV inclusion of the cubic term 
improves the quality of approximation, and at 1/r = 100 GeV also the term 0(a^J is needed. 
Let the truncution at n = n + result in the closest approximation of c% from above and 5V + {r) 
denote the difference between the exact and the approximate values of the potential Vri(r). If 
the contribution of the term for n = n + is not included, i.e. the series in (|29|) is truncated at 
n = n + — 1 the closest approximation from below is obtained. Let 5V_ (r) be the corresponding 
difference of the potentials. The values of the functions 5V±(r) are also given in Table [I]. One 
can estimate the theoretical uncertainty of the perturbative approach by considering the change 
in the value of the potential 

SV(r) = 5V + - 5V^ (31) 

corresponding to truncating the asymptotic formula at n = n + — 1 and n = n + . It is interesting^ 
that 5V(r) does not change drastically with r and is of the order of Ar. 

In perturbative calculations of the energy levels of a QQ system in this model, the truncation 
of the asymptotic series for the potential leads to an uncertainty in the perturbative pole mass 
rriQ of the heavy quark. This uncertainty is of the order of 5m ~ 5V(ro)/2, where ro denotes 
the typical size of the bound-states. This estimation is obtained assuming that the shift in the 
binding energy due to the change of the coupling a H is equal to — 5V(r ), which would be the 
case for the Coulomb potential. The shift in the binding energy is then compensated by shifting 
the masses of Q and Q by 5m. One may think that for a stable quark Q the uncertainty in 
mg can be reduced by comparing the QQ bound states of different sizes, thereby disentangling 
the correlation in the dependences on the coupling and mass. However, as it has already been 
mentioned the variation of 8V(r) with r is quite moderate. Moreover at larger distances non- 
perturbative effects become important. So, in fact it is difficult to reduce the uncertainty in 
the mass determination significantly.^ 

It may appear too pessimistic to estimate 5m from 5V. One can argue that instead of 
considering the truncated asymptotic series for ct H the exact expression (^fj) should be used in 
the Schrodinger equation for the QQ system. Then mQ can be determined from the mass of 
the QQ bound-states without any theoretical uncertainty. This argument is based, however, on 
the complete knowledge of the function a n (q) including the range of its argument where non- 
perturbative contributions are very important. In purely perturbative calculations 5m has to 

5 This can be shown easily for a small cur,: Since asymptotically /„ ~ 2™, it follows that n min ~ 27t/(/3oO!r) ~ 
m(l/(Af{r)). Then one can estimate the size of the last term included or rejected in eq. ( p9| ) to be Sa^ ~ Art, 
which leads to SV(r) ~ Ar. 

6 Behind this argument for the Richardson potential lies a following corresponding perspective in QCD. In 
principle one is able to determine the two independent fundamental parameters of QCD, tuq and a s , from 
a sufficient number of different physical observables Oi(rriQ, a s ) by exploiting their different dependences on 
these parameters. When we express the observables in perturbative expansions of a s , however, we encounter 
uncertainties of the order of Aqcd ~ 200 MeV originating from truncations of the asymptotic series. Then, if 
we try to extract the values of tuq and a s from these observables, the uncertainties should be attached to both 
of these paramters which cannot be reduced within a purely perturbative approach. 
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Figure 2: The momentum-space couplings cxy{q) for different values of a-jjg{m z )- 

be interpreted as a perturbative theoretical uncertainty in uiq due to the asymptotic character 
of the series (|29|). 

It seems reasonable to expect that the qualitative features of the asymptotic expansion (f20|) 
should be similar for the Richardson potential and for the QCD potential. This is a useful 
assumption because only limited information is available in the latter case. In perturbative 
QCD the coupling a v is expressed as an asymptotic series in a-^-g and only the first three 
terms are known P, |3[]. In Fig. a v (q) is shown for rif = 5 active flavours and q between 10 
and 100 GeV. The solid line corresponds to a-jjg{mz) = 0.118 and the dashed lines to changing 
this input value by ±0.003 which is the present error of this parameter ||. In QCD the first 
three coefficients of the j3 function (|7|) are also known JED), 0, [|: 

Pv,o = A)j 

38 

P Vtl = /3i = 102-—n f , 

2239 377 9 , / 3tt 2 \ , 
(3 v2 = i 85 4___ n/ + _ n/ 2 + 7r 2^ 18 _^j (33 _ 2n/) 

+ (726-^n / + i°V) C(3). (32) 
Using this information one can write a formula for a v (l/r) which is accurate up to fifth order 
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a. 



y {l + 2(3 2 f 2 ^ 2 + (5(3 p 1 f 2 + 6Po 3 f 3 ) a v 3 + 

(3 A 2 + 6 A/2) / 2 + 26 A 2 A is + 24 /5 4 / 4 ] a v 4 } + 0( 



(33) 



where a v = a v (q = 1/r') and a v = a v /(47r). Let us remark that the term quadratic in a v is 
absent because it has been absorbed by shifting the argument of the coupling a v . After simple 
algebra the following formula is derived: 
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1 - \ n -^ n f) y av + 
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7T 
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3 J 9 J 27 
The numerical values of the coefficients in the above expansions for 71/ = 5 are equal to: 



a. 



a v + 1.22454 a v 3 + 5.59618 a/ + 32.2015 a v , 5 + C(a 



V )l 



whereas for nj = 4 and n/ = 3 one obtains 



a v + 1.44676 a v 3 + 7.38182 a/ + 46.3717 a v 5 + 0{a x 

+ 64.0389 a v 5 + 0(a v 



a v + 1.68750 a v d + 9.45282 ct y L 



(34) 



(35) 



(36) 
(37) 



respectively. The numerical coefficients in the relations (|3q)-(|37|) are large. The asymptotic 
character of the expansion is evident and only for very short distances (1/r > 30 GeV) all five 
terms should be kept. At such distances the number of active flavours is nj = 5. The 0(a 4 ) 
and 0{aP) terms in ([35]) are equal for a v = 0.1738 and their contributions to 6t v are 0.0051. 
This implies that the truncation of the asymptotic series leads to an uncertainty of about 2.8% 
in a v . The value of a v corresponds to 1/r = 32 GeV, a distance which is not very different 
from those probed in ti production near the threshold. Phenomenological consequences of 
the above observation are discussed elsewhere ||11||. At the distances probed by bb states, the 



corresponding values of a v are so large that the quartic term in the expansion (|34|) must be 
rejected. Thus it would be simply inconsistent to keep other contributions of this order like 
those in the relation between a v and a-jjg even if this relation were known. Our analysis leads 
to a rather surprising conclusion that in the framework of purely perturbative QCD, theoretical 
uncertainties cannot be reduced below the present level ! (This is consistent with the discussion 
given in ||.) If we accept such a radical point of view, we can proceed even further by observing 
that the cubic term in (|34]) depends only on j3 and is the same for the QCD static potential 
and for the Richardson potential. The latter case suggests that it is not obvious at all if the 
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inclusion of the cubic term improves the accuracy of the formula for a v . It may well be that a 
more precise answer is obtained if the cubic term is also rejected. Without extra information 
beyond the purely perturbative approach, we cannot answer the question: to reject or not 
to reject. Clearly, keeping the cubic term implies a stronger attraction between the quarks. 
Thus, if m& is determined from the energy of the bb ground state, the value of mj, obtained 
using a perturbative calculation and the potential including the cubic term is larger than the 
value corresponding to rejecting this term. The difference in is of the order of 100 MeV 
and our analysis indicates that attempts to reduce this error to a much smaller value within 
perturbative QCD may be inconsistent. 
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